Thin Ohmic or superconducting strip with an apphed ac electric current 
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The complex impedance, currents, and electric and magnetic fields are calculated as functions 
of resistivity and frequency or London depth for a long thin strip with applied ac current. Both 
Ohmic and superconducting strips are considered. While the inductance per unit length of the 
strip depends on the strip length logarithmically, the sheet current, magnetic field, resistance, and 
magnetic susceptibility are independent of this length. It is found that the enhancement of resistance 
by the skin effect in thin Ohmic strips is much weaker (logarithmic) than in thick wires. 

PACS numbers: 74.25.Nf, 73.25.+i, 74.78.Bz, 74.78. Db 



The distributions of the electric and magnetic fields in- 
side and around a long wire carrying an applied ac cur- 
rent la is a standard problem in electrodynamics when 
the wire has circular cross section In this case, the so- 
lutions of Maxwell equations away from the current con- 
tacts depend only on the radial coordinate r and are eas- 
ily obtained for both Ohmic wires with resistivity p and 
superconducting wires in the Meissner state with Lon- 
don penetration depth A. However, when the wire has a 
non-circular cross section, e.g., for flat strips, the prob- 
lem becomes difficult. While the theory of thin strips in 
an applied ac field with no applied current is known since 
long time - it is simpler since the strip length drops out 
from all final results - apparently, no analytic or trans- 
parent numerical solutions are available for the even more 
important response of thin strips to applied ac current. 
This response is required, e.g., to understand SQUID sys- 
tems or to analyze the important experiments on tem- 
perature dependent magnetic field profiles of supercon- 
ducting strips In particular, the dependences of in- 
ductance, resistance, and field and current distributions 
on A, p, ac frequency, and strip length, are of principal 
interest but are still unknown. 

In this paper we settle this fundamental question and 
obtain these dependences for a current-carrying thin strip 
of width 2a, length I ^ a, and thickness d a, filling the 
volume < a, \y\ < d/2, \z\ < I. In this basic case, the 
magnetic field around the strip depends only on the sheet 
current J{x) = J^^^2 jX^j?/) "^2/1 which like the current 
density j and total current la flows along z. Interestingly, 
the inductance per unit length of a long strip depends on 
its length I although the field and current distributions 
in the strip and its resistance and magnetic susceptibility 
are independent of /. 

This strip example has the additional advantage that 
the nontrivial magnetic field distribution at the surface of 
the conductor can be measured by magneto-optics using 
flat indicator films or linear arrays of Hall probes [3, 0] ■ 
Note that with cylindrical wires there is no analog to such 
field profiles since the wire surface r = a exhibits constant 
and material-independent field magnitude Ia/2T:a. 

To solve the strip problem with minimum effort we de- 



rive an equation of motion for the sheet current J(x,t), 
which in general is driven by a time-dependent perpendic- 
ularly appfied magnetic field Ha{x,t)\\y and/or an elec- 
tric current Ia{t) applied at the ends of the strip, far 
away from the considered section. This method implic- 
itly accounts for the complicated magnetic field around 
the strip, which thus does not have to be cut off or ap- 
proximated as in other numerical methods. For strips 
with no applied current in a perpendicular Ha(t), inte- 
gral equations for J(x,t) and j{x,y,t) were obtained in 
Refs. |5| for thin and thick strips with linear or nonlinear 
resistivity, assuming zero London depth A = 0. But ac- 
tually, from the linear results of Refs. one can obtain 
the equations and results also for a superconductor strip 
in the Meissner state with arbitrary A, namely, by replac- 
ing the Ohmic (real) resistivity p by an imaginary value 
p{ijo) — iujpoX^ . The circular frequency lu then drops out 
from the final (static) results. 

General equations for thin and thick strips with both 
applied Ha and la and with finite A were derived in Ref . 
using the rather general voltage-current law valid inside 
superconductors with or without vortices. 



E^E,{v,J,li) + ^io\\v)^J/^t. 



(1) 



In it i?t,(r, j, B) = /9„(r,j, B)j' is the electric field caused, 
e.g., by moving vortices, and B is the magnetic induction. 
The second (London) term describes the acceleration of 
the massive charge carriers (Cooper pairs) by the electric 
field. The flux-flow term can be linear or nonlinear, 
e.g., pi, = pff = const'-B for free fiux fiow, or py oc 
for thermally activated depinning of vortices, where n 3> 
1 is the flux-creep exponent. 

For the particular case of a thin strip with applied 
current la but no applied field Ha (or with large dc 
background field), the sheet current has the symmetry 
J{—x) — J{x). The equation for J{x,t) then reads 

^^^^ = ti^'£dx'K{x,x') [Ea{t) - E4x\t)] . (2) 

Here E^{x',t) = Ey{x',j,B) with j = J {x' , t) / d{x') and 
B = poH{x',t) the magnetic induction in the strip. The 
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spatially constant electric field Ea{t) (along z, like j and 
J) formally drives the currents; it will drop out later 
when the solution is expressed in terms of /a(i)- Fi- 
nally, the integral kernel K{x, x') is the inverse of a kernel 
K~^{x^ x') and is defined by 

[ dx"K{x, x") K-\x", x') = S{x - x'), 



K-\x,x') = —\n 



2tt 



+ A{x) S{x - x') , (3) 



with A(x) — \^{x)/d{x) the effective magnetic penetra- 
tion depth of a superconductor film with thickness d < X. 
Note that Eqs. (2) and (3) apply also to inhomogeneous 
strips with A, p^, and d depending on sufficiently 
smoothly (over length scale ^ d\^ _This property can 
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be used to simulate edge currents |3, |J| caused by a geo- 
metrical barrier 0,0 for flux penetration. 

For a strip with Ohmic resistivity p and applied ac 
current Ia{t) = lao^^pi'i^t) one has A = and Ey = 
Pvj — {Pv/d)J, or one may formally put Ey = and A = 
—ip/{p,oUjd). Both methods yield the same equation for 
J{x,t) — Jq{x) exp{iu)t) or E{x,t) = pJ{x,t)/d driven 
by Ea{t) = Eaoe^p{iujt), 



Jo{x) 



iEaQ 



Jo 



dx' K{x,x'). 



(4) 



The integral kernel K{x, x'), Eq. (3), is now complex due 
to the imaginary A. Dividing this by the current ampli- 
tude lao = 2 Jq dx Jo (a;), the electric field amplitude Eao 
drops out: 



Jo{x 



. dx' K{x,x') /2 / dx dx'K{x,x'). (5) 

laO Jo I Jo Jo 

The magnetic field caused by this sheet current is 
li{x, y,t) = Ho(a:, y) exp(jwi) = {Hx,Hy) with 



Ho(a;, y) — Vx z dx 

J —a 



In- 



{x — x')^ + y 



(6) 



Interestingly, the profiles Jo (a;) and Ho(a;,y) obtained 
numerically from Eqs.(4)-(6) do not depend on the strip 
length I, though the kernel K{x,x') looks different for 
different ratios I /a. This may be understood by physical 
arguments, but to prove it mathematically appears to be 
difficult. See also the similar theory for a double strip 
0. However, the complex resistance or impedance Z = 
iujL + R per unit length does depend on I logarithmically 
{R = real resistance, L — real inductance). Noting that 
Eal is the voltage drop along the strip, one obtains 



Z 
J 



EaO 



iLupo / 2 1 dx J dx'K{x, x') . (7) 



We evaluate the integrals (2)-(7) using unit length a = 
1 and introducing a grid Xi {0 < Xi < 1, i — 1,2, . . . , N) 
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FIG. 1: The complex amplitude of the sheet current Jo{x) 
in superconducting and Ohmic strips with applied ac cur- 
rent, Eq. (9). Shown are the real (in-phase) and imaginary 
(out-of-phase) components in units of lao/a for various 2D 
penetration depths A or skin depth depths 5 (|A| = 5'^ /2d). 



with weights Wi such that for any sufficiently smooth 
function f{x) one has j^f(x)dx « Tl!i=i'^if{xi)- ^ 
grid with grid points densely spaced near the strip edge 
a; = I, is obtained by introducing an auxiliary grid 
Ui — [i — and then putting, e.g., Xi — {iui — m?)/2, 

Wi — 3(1 — uf)/2. This amounts to a substitution of 
variables in the integral, Jq f{x) dx = f[x{u)] x'{ii)du. 
The integral kernel K{x, x') now becomes a matrix 



K,. 



Qij — 



{wjQij + A{xi) Sij) S 



— In 

27r b.? - 



1 , 7r;2 
— In 

27r XiWi 



(8) 



where Sij equals 1 when i — j and otherwise. The diag- 
onal term Qa in Eq. (8) is chosen such that the numerical 
error decreases with a high power of the grid number N, 
say, as N^"^ or N~^, depending on the chosen grid, while 
any different choice gives a larger error oc 7V~^. The 
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FIG. 2: The complex amplitude of the perpendicular mag- 
netic field Hoy(x,Q) in the plane of superconducting and 
Ohmic strips with applied ac current, Eq. (10). Shown are 
the real and imaginary parts in units of lao/a, for various |A|. 



current and field profiles are then obtained as 

N , N N 



laO 
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k=l 1=1 



(9) 
(10) 



where V.^ means the numerical derivative d/dx; this ex- 
pression (10) is much more accurate than the direct com- 
putation of HQy{xi) from Jo(a;i) by Ampere's Law using 
a matrix with singular terms oc 1/ {xi — Xj). 

The profiles Jo{x) and Hoy{x,0) for various A are 
depicted in Figs. 1 and 2. For Ohmic and super- 
conducting strips one has the same limits [loj : at 
|A| ^ a (ideal screening) one has Jo{x) = (/ao/''r)(a^ — 
a;2)-i/2^ Hoy{x,0) = for \x\ < a and Hoy{x,0) = 
(/qo/Stt) sign(a;)(a;^ — a^)^-^/^ for |a;| > a, and at 
A| » a (weak screening) Jo{x) — (/qo/So) — const 
and Hoy{x,0) = (/ao/27r) sign(x) In |(a - \x\)/{a + \x\)\. 
The imaginary (out-of-phase) components of Jo{x) and 
Hoy{x,0) are zero in these limits, but are maximum for 
Ohmic strips at some x dependent value of |A| ~ S'^/2d 



where S is the skin depth. 

The complex impedance, Eq. (7), is calculated as 



N N 



Z = iu)L + R — iiOfiol / WkKki . 



(11) 



fc=i 1=1 



For superconducting strips this yields R 
inductance (partly obtained analytically) 
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(12) 



with a{0) = and a{oo) = 3/2-ln4 = 0.11371. To good 
approximation, the first term is the geometric inductance 
Lm and the last term the kinetic inductance Lk, while 
the small middle term cx a is shared by both and L^. , 
which slightly depend on the current distribution. Fig. 1. 
For Ohmic strips we obtain Z = iluL + R, 
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(13) 



with |A| = p/{^QUjd) = 5^ /2d and i?strip = pl/2ad. The 
functions a, (3 and 7 are shown in Fig. 3. One has 
/3(0) = (ideal screening), /3(oo) = a(oo) = 3/2 — In 4 
(uniform current), 7(00) = 1, and for |A| ^ a, 7 « 
1 + (2/7r2)ln(0.14a/|A|) (skin effect). The same 7 re- 
sults from Ohmic dissipation, 7 = 4a JJ" \ Jq[x) / Iao\^dx, 
Eq. (5). Good fits for all A are (see Fig. 3) 

a « a(oo) ( tanh[ 0.435 ln( A/0.074 a) ] -h l)/2, (14) 

13 « a(oo) V(tanh[0.921n(|A|/0.15a)] l)/2, (15) 
7 « 1 + (4/3^2) i^[i^(o.l4a/|A|)3/2]. (16) 

For comparison we give here also the correspond- 
ing expressions for a uniform round wire of radius a, 
with applied current Ia{t) = /ao exp(itt't), current den- 
sity j{r,t) = jo{r) exp{iujt), and magnetic field H = 
(fHo{r) exp{iujt), where ip is the unit vector in azimuthal 
direction. Outside the wire (r > a) one has simply 
Ho{r) = /ao/27rr, irrespective of the material properties 
and of uj. Inside this wire (r < a) one has 



Hoir) 
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IgQ hjr/X) 
27ra /i(a/A) 
laO Io{r/X} 
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djoir) 



dr ' 
1 5[ri/o(0] 



/i(a/A) 2A r dr 



(17) 
(18) 



where /o(a;) and Ii{x) are Bessel functions. When A is 
real, Eqs. (17) and (18) are well known London solutions, 
and with imaginary A^ — {p/iuj^lq), or complex A^^ = 
(1 -I- i)(w/xo/2p)^/^ = (1 -I- i)/5, they describe the skin 
effect in Ohmic wires with skin depth 5 = (2/9/a;/io)^^^- 

The complex resistance Z — icuL + R and self- 
inductance L of this cylindrical wire may be obtained 
from the energy balance . In terms of real la (t) and 
Ea(t), the input power laEJ is the sum of the tempo- 
ral change of the magnetic field energy outside the wire. 
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FIG. 3: The functions a, /3, and 7 of \A\/a entering the 
inductance and resistance of a strip, Eqs. (12) and (13) (dots). 
The solid lines show the approximations Eqs. (14)-(16). 



{d/dt)nol C"" 2TTr{Ia/2TTrfdr ^ {d/dt)LJl/2 defining 
the "external inductance" 



wires with A = 5/{l + i) making the second term com- 
plex. For superconducting wires with A ^ a this yields a 
kinetic inductance = hqX^I/tto? . For Ohmic wires 
with large 5 ^ a, expansion of (21) adds a constant 
j3 = 1/4 to h\{l/a) (from the "inner inductance") and 
a factor 7 = 1 + a'*/48(5'* « 1 to i?wiic- For smah skin 
depth (5 ^ a, this constant \s (3 — S/a and the factor is 
7 = i?/i?wiro = a./5 :s> 1, i.e., the resistance R is strongly 
enhanced by the skin effect. 

In conclusion, the dependences on A, p, uj, and I, of 
the inductance L and resistance R of superconducting 
and Ohmic strips with apphed ac current are obtained, 
Eqs. (12) to (16), and the current and magnetic field pro- 
files are depicted. As compared to cylindrical wires, in 
strips the enhancement of R by the skin effect is much 
weaker, R/Rstnp = 7 « H-0.2 \n{0.3ad/S^) for small skin 
depth S a and d < 6. Superconducting strips have a 
kinetic inductance Lk ~ /io^A/2a = /ioA^//2ad. The ge- 
ometric inductance of both Ohmic and superconducting 
strips is similar to that of cylindrical wires; it is domi- 
nated by the length-dependent factor ln(2Z/a) that origi- 
nates from the magnetic field energy outside the conduc- 
tor and does not depend on material parameters. This 
factor ln(2Z/a) results from Eq. (7) without having to 
introduce an outer cut-off here. The strip length / en- 
ters when the integral kernel K{x, x'), Eq. (3), is derived 
by integration of the 3D Laplacian kernel l/47r|r — r'| 
over z' , assuming z independence of J. But, of course, 
the exact value of the cut-off length in the logarithm de- 
pends also on the return path of the current. 

I thank J. R. Clem, G. P. Mikitik, and E. Zeldov for 
helpful discussions. This work was supported by the Ger- 
man Israeh Research Grant (GIF) No G-705-50.14/01. 
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with outer cut-off radius rout ~ I, and the power dissi- 
pated or stored inside the wire (in the magnetic field and 
kinetic energy of the Cooper pairs). Using the Point- 
ing vector, this power may be expressed in terms of the 
surface values = E{a,t) and = H{a,t), thus 
laEJ = {d/dt)LJl/2 + HsEs ■ 2Tral. Inserting here 
Hs = Ia/2TTa and dividing by la we obtain 



Ea{t)l^ Lja{t)+Es{t)l. 



(20) 



Writing the linear Eq. (20) in complex notation with 
Eao = Zlao, Es = pjo(a)exp(iu;i), jo{a) from Eq. (18), 
and i?wiio = pl/T^a^, we find 

^ = ^"27l"a+^-'-2A7rRA)- ^2^) 

Formula (21) applies both to superconducting wires with 
real London A (and i?wiro = iuj^QX^l/ira'^) and to Ohmic 
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